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Abstract. We present data on long-range multiplicity correlations in ultra- 
relativistic heavy ion collisions at the top RHIC energy {^/sWn — 200 GeV) 
. from the STAR experiment. The data shows a long-range multiplicity corre- 

lation extending across a gap of 1.6 units in pseudorapidity. The data is over 
. predicted by a multiparticle production model with hadronization of indepen- 

dent strings, or fusion of two color strings. This can be interpreted in terms of 
additional dynamical reduction in the number of particle sources. 
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pH I 1. Introduction 

The measurement of particle correlations has been suggested as a method to search 
, for the existence of a phase transition in ultra-relativistic heavy ion collisions [ 

' Q Ql 0] ■ If the quark-gluon plasma (QGP) is formed in these collisions, the existence 

or absence of particle correlations could lead to a determination of the presence of 
partonic degrees of freedom. Multiparticle production at low pt in 200 GeV is well 
described by string hadronization models, such as the Parton String Model (PSM), 
a Monte Carlo realization of the Dual Parton Model (DPM), or Quark-Gluon String 
Model (QGSM), which includes fusion of overlapping color strings [|2 HI Ej. The 
color strings are formed between the projectile and target nuclei. These strings 
cover a small area in the transverse collision plane that contains the color field 
created by colliding partons. The number of color strings is directly related to the 
atomic number or energy of the colliding nuclei increases, so does the number of 
color strings. The strings will then begin to overlap in the transverse plane, forming 
clusters. The PSM (an evolution of the String Fusion Model [I3), incorporates the 
fusion of two strings produced from soft collisions. 
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The transverse area of the strings goes as -; — ttt, so there is a low probability of 

fusing hard strings (strings with large mean transverse momentum, (pt))- The fused 
string conserves energy-momentum and color charge. Two immediate implications 
of string fusion are a reduction in overall particle multiplicity and an increase in 
(pt) of produced hadrons. Both these effects are observed in Au+Au data at RHIC 
[IH]. Hadronization of the strings proceeds with a Schwinger probability dependence 
[ EI ■ There is no a priori assumption of a phase transition in either the DPM or 
PSM. However, the clustering of strings may eventually lead to the formation of 
a macroscopic cluster at a certain critical density. This would indicate the onset 
of the percolation phase transition [1^1 E|- In this context, string fusion can be 
interpreted as an intermediate stage leading toward the formation of deconfined 
quark-gluon matter 

A linear relationship has been found in high-energy colliding hadron experi- 
ments between the multiplicity in a forward rj region (Nf) and average multiplicity 
in a backward rj region fA^h')[ll2): 

< NbiNf) >^ a + bNf (1) 

The coefficient b is referred to as the correlation coefficient and is a function of 
incident energy and atomic number. It can be expressed in terms of the expectation 
value [E]: 

_ < NfN, >-<Nf><m> _ 

<N]>- <Nf>^ ~ D)f 

where D^j and are the backward-forward and forward-forward dispersions, 
respectively. This result is exact and model independent 

Short and long-range (in rapidity) multiplicity correlations are predicted as a 
signature of string fusion [^J When strings fuse, a reduction in the long-range 
forward-backward correlation is expected. The existence of long-range multiplicity 
correlations may indicate the presence of multiple partonic inelastic collisions. These 
correlations arise from the superposition of a fluctuating number of strings, such 
that [ inj : 

Dlf =< TV/iVfo > - < Nf >< Nb >cx [(< 71^ > - < n >^)] < Nqj >< iVoh > 

(3) 

with {{p?) — {n) ) the fluctuation in the number of inelastic collisions and (A'o/) , (A^ofo) 
are the average multiplicity produced from a single inelastic collision. Therefore, 

should be sensitive to the presence of long-range multiplicity correlations. 

In this paper we discuss the results from 200 GeV Au-|-Au collisions for all 
charged particles withpr in the range from 0.1 to 1.2 GeV/c, to ensure a sampling of 
soft particles only. To eliminate short-range correlations, a gap in pseudorapidity {rj) 
of 1.6 units is considered [|^. The forward pseudorapidity interval was 0.8 < rj < 1.0 
and the backward was — 1.0 < < —0.8. 
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2. Data Analysis 

The data utilized for this analysis is from year 2001 (Run II) ^snn — 200 GeV 
Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC), as measured by the 
STAR (Solenoidal Tracker at RHIC) experiment The main tracking detector 

at STAR is the Time Projection Chamber (TPC) [EI. The TPC is located inside 
a solenoidal magnet generating a constant, longitudinal magnetic field. For this 
analysis, data was acquired at the maximum field strength of 0.5 T. All charged 
particles in the TPC pseudorapidity range 0.8 < \r]\ < 1.0 and in the pr range 
0.1-1.2 GeV/c were considered. The range was chosen to eliminate the effects of 
short-range correlations from sources such as cluster formation, jets, or resonance 
decay. The collision events were part of the minimum bias dataset. The minimum 
bias collision centrality was determined by an ofiline cut on the TPC charged particle 
multiplicity within the range 0.5 < \r]\ < 0.5. The centralities used in this analysis 
account for 0-10, 10-20, 20-30, 30-40, 40-50, 50-60, 60-70, and 70-80% of the total 
hadronic cross section. An additional offiine cut on the longitudinal position of the 
collision vertex (vz) restricted it to within ±30 cm from z = (center of the TPC). 
A combined two million minimum bias events satisfied these requirements and were 
used for this analysis. Corrections for detector geometric acceptance and tracking 
efficiency were carried out using a Monte Carlo event generator and propagating 
the simulated particles through a GEANT representation of the STAR detector 
geometry. 

In order to eliminate the effect of statistical impact parameter (centrality) 
fluctuations on this measurement, each relevant quantity {{Nf), (Nj,), (Nf) , and 
(NfNi,)) was obtained on an event-by-event basis as a function of STAR reference 
multiplicity, N^h- A linear fit to (Nf) and (iVf,), or a second order polynomial fit 
to (iV/)^ and {NfNb) over the range < Nch < 600 was used to extract these 
quantities as functions of Nch- Tracking efficiency and acceptance corrections were 
applied to each event. These were then used to calculate the backward- forward and 
forward- forward dispersions, D^j and D^f, binned according to the STAR centrality 
definitions and normalized by the total number of events in each bin. 

3. Results 

Fig. n presents the mean multiplicity {Nf), (iVfc), in the forward and backward 
pseudorapidity region (0.8 < |ry| < 1.0) along with (iV/)^ and (NfNb) as functions 
of reference multiplicity, Nch- The results are calculated event-by-event for every 
unit Nch- The {Nf) and (Nt) demonstrate a linear dependence over much of the 
centrality range. Only statistical uncertainties are present and are less than 1%, 
except for the highest multiplicities. These results are uncorrected for tracking 
efficiency and acceptance. Due to statistical limitations, it is not possible to apply 
corrections for every value of Nch- One value for the correction, calculated for 
each centrality, is applied to the resultant values {N f) , (Nb) , {N f)'^ , and {NfN},) 
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Fig. 1. {Nf),{Nb), {Nf)", and {NfNb) as functions of STAR charged particle 
reference multiplicity (Nch)- The quantities are found for event-by-event for all 
values of Nch- 



for every event in that centrality. Therefore, all events falling within a particular 
centrality have the same correction. 

Fig. 121 shows the results for D^f, and the correlation coefficient, b, as 

a function of Nch for the 8 centrality bins. The presence of long-range multiplic- 
ity correlations are evident from D^j and b. The growth of D^j as a function of 
Nch is consistent with an increasing long-range correlation from peripheral to cen- 
tral heavy-ion collisions, corresponding to a greater number of fluctuating strings. 
The linear evolution of D^j as a function of Nch represents a (short-range) auto- 
correlation of particles within the same window. The b coefficient represents the 
long-range correlation normalized to the short-range correlation in the forward 77 
window. Were there no long-range correlations [D'^f = 0), then b would vanish. 
The vertical error bars in Fig. El are statistical and systematic uncertainties. The 
statistical errors are smaller than the data points. The systematic errors are deter- 
mined by varying cuts on \vz\, number of fit points on track in the TPC, fit ranges 
for (TV/) , (Nb) , {Nf) , and {NfNb), etc. Horizontal error bars are the RMS value 
from the mean multiplicity for each centrality. 
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Fig. 2. (a) D^j, (b) Djp and (c) b from Au+Au collisions at ^/snn — 200 
GeV, as a function of the STAR reference multiplicity, Nch- Vertical error bars 
represent statistical and systematic uncertainties, horizontal, the RMS of the mean 
multiplicity distribution in each centrality bin. 



The comparison of and as calculated from the Au+Au data to that 
from the Parton String Model (PSM) with the string fusion on or off is presented in 
Figs. 121 and 01 respectively. The PSM (with two string fusion) minbias multiplicity 
distribution closely matches that of the corrected Au+Au data. It also describes 
the (pt) enhancement, particle ratios, strangeness production, etc., seen in Au+Au 
collisions at RHIC [El- As such, the centrality cuts in the PSM correspond to the 
same percentage of the minimum bias cross section as in the data and are plotted as 
a function of the number of participant nucleons in the collision (Npart), calculated 
from Monte Carlo Glauber model [El- The statistical and systematic uncertainties 
in Figs. Oland^are the same as those described for Fig. 

Figs. |Sland0]show good agreement in peripheral collisions between data and the 
independent (no string fusion) or collective (with string fusion) model. In central 
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and (b) Dj:j as a function of the Npart in 200 GeV Au+Au collisions, 



Fig. 3. 

compared to the PSM with no string fusion (independent strings). There is good 
agreement for both quantities in peripheral collisions. A large discrepancy exists in 



central collisions, suggesting an underlying dynamical mechanism that reduces 



and Djj. 



Au+Au collisions, 1?^^ with the independent string description deviates from the 
data, but shows excellent agreement for fusion of two soft strings. This confirms the 
agreement in multiplicity between the PSM (with two string fusion) and data, as 
D'jj: is a multiplicity effect. However, there is a large discrepancy in D^^ for central 
collisions for both the independent and collective PSM, compared to the data. This 
discrepancy is greater for the case of independent strings (Fig. ISJ. This suggests an 
additional, dynamical reduction in the number of particle sources in central Au+Au 
collisions, greater than that provided by the fusion of two soft strings (Fig. 

4. Conclusions 

We have presented data on long-range multiplicity correlations of soft particles in 
200 GeV Au+Au collisions. A long-range multiplicity correlation is seen for all cen- 
tralities and increases from peripheral to central collisions. Comparison to the Par- 
ton String Model (PSM) using both interacting and non-interacting strings shows 
good agreement for peripheral collisions, but large differences in central collisions. 
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Fig. 4. (a) and (b) Djj: as a function of the Npart in 200 GeV Au+Au 
collisions, compared to the PSM with string fusion (collective strings). There is 
good agreement for both quantities in peripheral collisions. The discrepancy in 
central collisions is reduced compared to Fig. |31 but still persists. This remaining 
difference hints that an additional reduction, beyond that offered by two soft string 
fusion, is required. 



The reduction in the long-range multiplicity correlation, D^j: , even compared to the 
collective PSM (fusion of two soft strings) , is a possible indication of an additional, 
collective effect in central Au+Au collisions at 200 GeV. 
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